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1. Introduction 

Lattice QCD has by now entered the era in which high precision calculations of physical 
quantities at the physical point including dynamical fermionic degrees of freedom become feasible 
(see for example the quantities summarised in the FLAG review [1]). This progress has been 
possible lai^gely due to substantial algorithmic improvements and the associated acceleration of 
the simulations allowing to simulate light quarks. Despite these successes there are quantities 
where a satisfactory control over systematics has not been achieved so far, among them quantities 
related to the structure of hadrons [2, 3]. For the pion electromagnetic form factor the dominant 
systematic uncertainty is due to the chiral extrapolation. The form factor has been worked out in 
chiral perturbation theory [4, 5] ixPT) to next-to-next-to leading order [6, 7] (NNLO) which can 
be used to guide the chiral extrapolation. However, the pion mass range for which ;fPT to a given 
order is applicable is not known a priori. 

In this proceedings article we present the status of our calculation of the pion form factor 
fnn{Q^) in the region of small space-like momentum transfers = —q- in Lattice QCD with two 
dynamical quark flavours. In this regime the form factor encodes information about the electro- 
magnetic charge radius (r^) of the pion. 



Determinations of the charge radius usually suffer from an intrinsic model dependence since the 
slope of the form factor is obtained from fits over a large Q^-interval assuming a particular model 
for the 2^-dependence such as pole dominance. Here we use suitably tuned partially twisted bound- 
ary conditions [8, 9, 10, 11, 12] to compute data points in the immediate vicinity of 2' = which 
enables the model-independent extraction of the charge radius. The data is then compared to xVY 
which is used to extrapolate the charge radius to the physical point. The continuum extrapolation 
is done by introducing terms that model the cutoff dependence into the chiral extrapolation. The 
systematic uncertainty which is introduced by relying on xVY for our range of pion masses is taken 
into account by checking the results from xVY with results from a polynomial extrapolation. Parts 
of the results were already reported in [13, 14, 15, 16]. 

2. Extraction of the form factor 

The electromagnetic form factor in two-flavour QCD is defined by 



In lattice QCD the matrix element can be extracted from suitable ratios of 2- and 3-point functions, 
in particular we use the ratio [12] 




(1.1) 




(2.1) 




(2.2) 
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i8 


ro/a 


lattice 


nin [MeV] 




Labels Statistic 


5.20 


6.15 (6) 


64 X 32^ 


455,351,299 


6.0-4.0 


A3-A5 ^(100) 


5.30 


7.26 (7) 


64 X 32^ 


552,412 


6.2, 4.7 


E4,E5 ^(100) 






96 X 48^ 


295, 254 


5.0, 4.2 


F6, F7 ^(250) 


5.50 


10.00(11) 


96 X 48"* 


625, 534, 424 


7.7-5.3 


N3-N5 ^(100) 



Table 1: Compilation of simulation parameters. 



where 

C2(r/2,o) 

2C3(f,0,0) 

ensures the physical normalisation fj^ji (0) = 1 and 

C2{t,p) = |2y>P--^-)(P(^x)P(0,x,)), 

.7, ,x 

Here P{x) and V^(x) are the pseudoscalar density and the (a) -improved vector current, respec- 
tively, and E{p) is the energy of the pion with momentum p. Up to corrections from excited states 
the ratio in eq. (2.2) is constant in Euclidean time and gives the desired matrix element of the vector 
cun^ent from eq. (2.1). In practice, the factor is obtained prior to the extraction of the matrix 
element via eq. (2.3). 

3. Simulation setup 

The calculations are based on the CLS ensembles that have been generated with two degenerate 
flavours of non-perturbatively ^(a)-improved Wilson fermions using the deflation accelerated DD- 
HMC [17] and MP- KMC [18] algorithms. The bare parameters and some of the basic properties of 
the lattices are given in table 1. We express dimensionful quantities in units of the Sommer scale 
ro [19], which was recently determined on the CLS ensembles [20, 21] with a physical value of 
ro = 0.503(10) fm [21]. 

All 2- and 3-point functions were computed using two hits of stochastic Z2 x Z2 wall sources [22, 
23, 24, 25] on two different timeslices separated in time by T /2. The only exception is ensemble 
A5 where we have used four source positions separated by r/4. We ai^e primarily interested in the 
region of small momentum transfers which is made accessible by partially twisted boundary condi- 
tions at vanishing Fourier momentum. For each ensemble we have used five different twist angles 
in the .x-direction, that are tuned to achieve a high resolution in the low region. In addition 
to the form factor we have also computed the renormalised PCAC quark mass m and pion decay 
constant F„. All quantities are fully i^(a)-improved, and we adopt the renormalisation constants 
from [26, 27, 28, 29, 21]. The results are listed in table 2. Note that our ensembles cover a large 
range of quark masses and lattice spacings, always satisfying nijiL > 4. 

Statistical errors have been estimated using the bootstrap procedure with 1000 bins. At present 
all fits for the extraction of the masses and matrix elements are uncon^elated. 
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Figure 1: Results for the form factor on all ensembles. The inset shows a zoom in the region of low 
momentum transfer. 

4. Results for form factor and charge radius 

The full set of results for the form factor is shown in figure 1. The plot demonstrates that 
partially twisted boundary conditions are a powerful tool to explore the small momentum transfer 
region of form factors. In particular, the high density of points in the low region enables the 
accurate and model independent extraction of the charge radius, as discussed below. 

The data is still subject to effects stemming from the finite box size, even though they are 
expected to be small, since nij^L > 4. We estimate their size in ;tPT. For the pion mass and the 
decay constant the relevant expressions have been derived in [30]. For the form factor with partially 
twisted boundary conditions, finite size effects are only known for the case where one of the twists 
is zero [31], i.e. either the initial twist 0,- = or the final twist Of = 0, and in the Breit frame [32] 
where 0, = —Of - As input for the ;^PT formulae NLO low energy constants (LECs) are needed. To 
become independent of external input we apply an iterative procedure: First we estimate the finite 
volume effects using the LECs from [33] as input. We then perform the extraction of the charge 
radius and the form factor and fit the data to NNLO ;^PT. This yields new LECs that are used for 
the next iteration. After two such iteration this procedure converged and we did not see any residual 
deviations in the LECs. 

From now on we constrain the data to those kinematical situations where finite size effects are 
known in ;^PT. To extract the charge radius we compai^e the results from linear fits to the results 
for the slope at = obtained from polynomial fits to order two and three up to some maximal 
value of Q^. The results from those different fits for ensemble F6 ai^e shown in figure 2 (left). We 
also show the results obtained from the vector pole dominance model, i.e. of fits to a single pole 
form, which has predominantly been used to extract the chai^ge radius from the data of the form 
factor. Surprisingly, the linear fit starts to deviate from the other functional forms at relatively small 
values of Q^, which shows that higher order terms become important already around (Qro)^ 0.1. 
Since we like to extract the charge radius using a fit to the same interval for all ensembles 
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ensemble 




roF^ 




yeff 


A3 

A A 

A4 
A5 


1.161(12) 

A onc/1 1 \ 
U.oyj(ll) 

0.761(11) 


0.280( 8) 

U.z4/(lz) 
0.251(14) 


0.090(3) 
0.0jz(3) 
0.040(2) 


0.73228( 7) 
0. /zooD(lz) 
0.72731(10) 


E4 

J— /»J 

F6 
F7 


1.406(16) 
1 04Rn3t 
0.752( 8) 
0.646( 7) 


0.287(10) 

?7in It 

0.254( 8) 
0.237( 8) 


0.128(5) 
078C4t 
0.041(2) 
0.029(1) 


0.74962( 8) 

0.741 19( 4) 
0.74030( 5) 


N3 
N4 
N5 


1.593(18) 
1.360(16) 
1.080(13) 


0.329( 7) 
0.304( 9) 
0.291( 8) 


0.188(5) 
0.139(4) 
0.091(3) 


0.77162( 3) 
0.76855( 3) 
0.76543( 3) 



Table 2: Results for basic quantities. 
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Figure 2: Left: The squared charge radius versus the maximal value of entering the fit at a pion mass of 
295 MeV (ensemble F6). Shown are the results from a linear fit, polynomial fits of degrees two and three, 
and the results of a fit to the vector pole dominance model (VPD). Right: Compilation of results for the 
charge radius in dynamical lattice QCD [35, 36, 37] and from experiment [34]. 



we have to use the fit to a polynomial of degree two with a maximal value of {Qr^f- = 0.22, as 
indicated by the vertical dashed line in figure 2 (left). In this case all ensembles have more than 
three data points for the form factor below this 2^-cut. Note, that this does not affect the model 
independent extraction of the charge radius, since we can explicitly check for the other ensembles 
that the results are consistent with the ones from a linear- fit to the very small 2^-region. The 
results for the charge radius are shown in comparison to results from experiment [34] and other 
collaborations [35, 36, 37] in figure 2 (right). 



5. Chiral and continuum extrapolation 

The remaining task is the extrapolation to the physical point. The natural framework for the 
chiral extrapolation is SU (2) xVY. Since we have found that to next-to-leading order (NLO) 
does not work for our range of pion masses [15] (for more details see also our upcoming pub- 
lication), it is necessary to use j^T to NNLO. The expression for the form factor at NNLO in 
continuum jVT has been worked out in [6]. 

For the continuum extrapolation we have to keep track of lattice artefacts, entering at order a^. 
To take these into account we include a term proportional to a^, which respects the normalisation 



5 



The electromagnetic form factor of the pion in two-flavour lattice QCD 



Bastian B. Brandt 



fnni'^) = 1- The resulting formulae then contain 6 free parameters and not all of them can be 
constrained sufficiently by a fit to the data for the form factor alone. To overcome this problem 
we perform a simultaneous chiral extrapolation of fmi{Q^) (or directly {r^)) together with F„ and 
ruji with respect to the isospin-averaged light quark mass m. An account of the fit formulae will be 
given in our future publication. ' The inclusion of the additional quantities further constrains some 
of the LECs that are less constrained by the nin and 2^-dependence of the form factor. In total the 
formulae including terms modeling cutoff effects now contain 14 free parameters and further input 
is needed to be able to obtain reliable results. We therefore decided to fix the LECs = —0.4(5) 
and I2 = 4.3(1) by TiTT-scattering [33]. We have verified explicitly that our fit results are not very 
sensitive to l\ and I2 by checking that the results for (r|) at the physical point do not change 
significantly even if we vary them by a factor two. 

To eventually perform the chiral extrapolation we have used three different procedures, i.e. 

• simultaneous fits to ;^PT to NNLO including fnniO^), Fn and m^j, 

• simultaneous fits to ;^PT to NNLO including {r^), Fji and m^r, 

• polynomial fits to (r|). 

The last two procedures make use of the model independent extraction of the charge radius dis- 
cussed in the previous section while the first procedure assumes that ;t;PT also provides the correct 
2^-dependence of the form factor. Thus the last two procedures are the preferred choices for ob- 
taining a model independent result for the charge radius at the physical point. The polynomial fits 
provide estimates for the systematic effects stemming from the assumption that ;^PT to NNLO pro- 
vides the correct functional form for the chiral extrapolation for the present range of pion masses. 
An estimate for the associated error can be obtained from the spread of the results of the three 
different procedures. 

Figure 3 shows an example for the chiral extrapolation using ;^PT to NNLO for {r\^, Fji and 
niji including only ensembles with a pion mass below 600 MeV. The fit describes the data well 
and yields an uncorrelated ;t^^/dof =1.3 with 11 free parameters. The result in the continuum (red 
band) and at the physical point for the charge radius is larger than the experimental value, but not 
significantly. The result for the the pion decay constant is consistent with the experimental result. 
The same is true for the LECs that are in good agreement with the typical results found in the 
literature. A more comprehensive discussion of the chiral extrapolations will be given in our future 
publication. 

6. Conclusions 

In this proceedings article we have given an account of our comprehensive study of the elec- 
tromagnetic form factor of the pion at small momentum transfers in lattice QCD. The high density 
of accurate data points close to 2^ = allows for the model independent extraction of the charge 
radius at finite lattice spacings. Finite volume effects have explicitly been studied and corrected 
using ;^PT to NLO and all quantities are fully ^(a)-improved. While ^PT to NLO fails to describe 

'The part of the formulae is identical with the ones given in [36]. 
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Figure 3: Results of a simultaneous fit to the data for ('"^), as extracted in the previous section, F-,^ and mj^ 
using x^T to NNLO with pion masses below {m-,i ro)^ — 2. The coloured dashed lines are the central values 
for the different lattice spacings and the red band and the red dashed Une represent the error band and the 
central value in the continuum. 



the data consistently in this range of pion masses, ;^PT to NNLO provides a good description. 
However, stabilising the fits is challenging, which agrees with what has been found in [35, 36]. 
Final numbers and detailed discussions will be given in our upcoming publication. 
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